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Boundary Harnack principle and gradient estimates for 
fractional Laplacian perturbed by non-local operators 

Zhen-Qing Chen* * Yan-Xia Ren^ and Ting Yang^§ 


Abstract 

Suppose d > 2 and 0 < /3 < a < 2. We consider the non-local operator = A“/^ -|- 5^, 
where 

S^f{x) := Inn A(d, ^ (/(^ + ^) “ 

Here b(x, z) is a bounded measurable function on x that is symmetric in z, and A{d, —fd) 
is a normalizing constant so that when b{x, z) = 1, becomes the fractional Laplacian := 

— (—A)^/^. In other words, 

C^f{x) := liin A{d, -jd) [ {f{x + z)- f{x)) fix, z) dz, 

^^0 J\z\>e 

where j^{x,z) := A{d,—a)\z\~^'^^^'> + A{d,—I3)b{x,z)\z\~^'^^^y It is recently established in 
Chen and Wang [12] that, when j^{x,z) > 0 on x there is a conservative Feller process 
having as its infinitesimal generator. In this paper we establish, under certain conditions 
on 6, a uniform boundary Harnack principle for harmonic functions of (or equivalently, of 
C^) in any K-fat open set. We further establish uniform gradient estimates for non-negative 
harmonic functions of X^ in open sets. 
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1 Introduction 

Let d>2, 0</3<a<2, and b{x, z) be a bounded measurable function on x with 
b{x, z) = b{x, —z) for x,z € M'^. Consider the non-local operator -|- S^, where 

S^f{x) := lim A{d, -/3) f (/(x + z) - f{x)) dz. 
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Here A{d,—j5) is a normalizing constant so that when b{x,z) = 1, 5^ becomes the fractional 
Laplacian := —(—in other words, A{d,—f3) = + /3)/2)/r(l — /3j2). 

Thus can be expressed as 


= lini [ {f{x + z)- f{x)) j\x, z) dz 


where 

■bf N A{d,-a) A{d,-I3)b{x,z) 

1 (“^.2) = -ppT^ +--■ 

Note that since b{x,z) is symmetric in z, for / G 

= f {f{x + z)- f{x) - V/(x) • 2 ;l{|^|<i}) j^{x, z) dz. 


Recently, the fractional Laplacian perturbed by a lower order non-local operator 5^, and its 
fundamental solution have been studied in Chen and Wang [12]. It is established there that if for 
every x G M'^, j^(x, z) > 0 (that is, b{x, z) > — for a.e. z G M'^, then has a unique 

jointly continuous fundamental solution p^{t,x,y), which uniquely determines a conservative Feller 
process on the canonical Skorokhod space ©([O,+oo),such that 


Ex 



f{y)p’'{t,x,y)dy, 


X G 


for every bounded measurable function / on The Feller process is typically non-symmetric 
and it has a Levy system {J^{x,y)dy,t) (see [T2l Theorem 1.3]), where 


j\x,y) := j^{x,y - x). 


( 1 . 1 ) 


When b takes constant value e > 0, has the same distribution as the Levy process Y + 
where Y and Z are rotationally symmetric a- and /3-stable processes on that are independent 
of each other. Moreover, two-sided heat kernel estimates for have been obtained in Chen and 
Wang m, while two-sided Dirichlet kernel estimates in open sets have recently been obtained 
in Chen and Yang m- We emphasize that the lower order perturbations considered in mills] 
and in this paper are of high intensity, in the sense that j^^^dz may be identically infinite. 
This is the case, for example, when b{x, z) is bounded away from 0. Such high intensity non-local 
perturbation is in stark contrast to the Feynman-Kac type non-local perturbation considered in 
13 [g ttn la [T6] . where the kernels used in the non-local perturbation are typically integrable. 
Under high intensity non-local perturbation, the potential theory for the perturbed generator 
can be different from that of A"/^. For example, the heat kernel of may not be comparable to 
that of A“/^; see [12] fo'^ details. 

In this paper, we investigate boundary Harnack principle and gradient estimates for non¬ 
negative harmonic functions of in open sets. Boundary Harnack principle (BHP) asserts that 
non-negative harmonic functions that vanish in an exterior part of a neighborhood at the boundary 
decay at the same rate. It is an important property in analysis and in probability theory on har¬ 
monic functions. We refer the reader to the introduction of mm for a brief account on the history 
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of BHP that started with Brownian motion and then extended to subordinate Brownian motions 
and to certain pure jump strong Markov processes. Since is typically state-dependent and its 
dual operator may not be Markovian, the BHP results in mm are not applicable to harmonic 
functions of In this paper, we establish uniform boundary Harnack principle on K-fat open sets 
for non-negative harmonic functions of by estimating Poisson kernels of in small balls. 

Gradient estimates for harmonic functions of elliptic operators and on manifolds have been 
studied extensively in literature, including the celebrated Li-Yau inequality. See m and the ref¬ 
erences therein and for a coupling argument. See also [21] for gradient estimates of the transition 
semigroups of some Levy processes using coupling method. Gradient estimates for harmonic func¬ 
tions for non-local operators are quite recently. In [3], a gradient estimate for harmonic functions of 
symmetric stable processes is obtained. Gradient estimates for harmonic functions of mixed stable 
processes were derived in [22|. It has recently been extended to a class of isotropic unimodal Levy 
process in |2n|. For gradient estimate for harmonic functions of the Schrodinger operator + q, 
see [1] for a G (1,2) and [T9| for a G (0,1]. The second main result of this paper is to establish 
gradient estimates for positive harmonic functions of C)'. As far as we know, this is the first gradient 
estimate result for non-Levy non-local operators. 

We now describe our main results in details. In this paper, we use as a way of definition. 
For a, 6 G M, a A 6 := min{a, b} and a\/ b := max{a, b}. Let jx — y\ denote the Euclidean distance 
between x and y, and B{x,r) the open ball centered at x with radius r > 0. For any two positive 

C 

functions / and g, f ^ g means that there is a positive constant c such that f < eg on their 
common domain of definition, and f ^ g means that c~^g < f < eg. We also write and “x” 
if c is unimportant or understood, li D C. is an open set, for every x,y ^ D, define 

6d{x) := dist{x,dD) and rD{x,y) := 6 d{x) + Sniy) + \x - y\. (1.2) 


It is easy to see that 

roix, y) X 5d{x) -f |x - y| x dniy) + \x-y\. (1.3) 

Denote by := inf{f > 0 : 0 D}, the exit time from D by . When there is no danger of 

confusion, we will drop the superscript b and simply write td for r^. 

Definition 1.1. A function / defined on is said to be harmonic in an open set D with respect 
to X^ if it has the mean-value property: for every bounded open set U d D with U C D, 


fix) = E, 



for x € U. 


It is said to be regular harmonic in D if (11.41) holds for U = D. 


(1.4) 


Denote by 9 a cemetery point that is added to D as an isolated point. We use the convention 
that := d and any function / is extended to the cemetery point d by setting f{d) = 0. So 
Ea; [/(A)1^)] should be understood as E^, [/(A.(!^) : To < -|-oo]. In Definition ll.il we always assume 
implicitly that the expectation in (11.41) is absolutely convergent. 
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Assumption 1 Suppose Mi, M 2 > 1. b{x,z) is a bounded function on x satisfying 

\\b\\oo < Ml and b(x, z) = b{x,—z) forx, 2 ;GM'^, (1.5) 

and there exists a positive constant Sq ^ [0,1] such that for every x,y € 

M^^r°{x,y) < J’^{x,y) < M2r°{x,y). ( 1 . 6 ) 

Here J^°{x,y) = A{d, —a)\x — + eQA{d, —(3)\x — y\~'^~^. Since J^°{x,y) depends only on 

\x — y\, we also write J^°(|x — y\) for J^°{x,y). 

Definition 1.2. Let K G (0,1). An open set D C is said to be K-fat if for every z G dD and 
r G (0,1], there is some point x € D so that B{x,Kr) C D H B{z,r). 

The following is the first main result of this paper. 

Theorem 1.3 (Uniform boundary Harnack inequality). Suppose Assumption 1 holds and D is 
a K-fat open set in with k G (0,1). There exist constants n = ri{d,a, (5, Mi) G (0,1] and 
Cl = Ci{d,a,l3,K,Mi,M2) > 1 such that for every zq G dD and r G (0,ri/2], and all non¬ 
negative functions u, v that are regular harmonic in Dn B{zq, 2r) with respect to and vanish in 
D'^ n B{zo, 2r), we have 

^<CM for,,y^DnB(.„r). 

v(x) v{y) 

We call the above property uniform boundary Harnack principle because the constants ri and 
Cl in the above theorem are independent of Eq S [0,1] appeared in condition (II.6p . We next study 
the gradient estimates for non-negative harmonic functions in open sets. We write or di for ^ 
and V for {dxi , • ■ • , dxd ) • 

Theorem 1.4. Let D be an arbitrary open set in Under Assumption 1, there is a constant 
C 2 = C 2 {d,a, P, Ml, M 2 ) > 0 such that for any non-negative function f in which is harmonic 
in D with respect to , V/(x) exists for every x G D, and we have 

\Vf{x)\<C 2 ^ . forxeD. 

1 A dD[x) 

For x = {xi ,■■ ■ ,Xd) G and 1 < i < d, we write x® for (xi, • • • , Xi-i, Xj+i, • • • , Xd) G . 
Assumption 2. Suppose there is i G {1, • • • , d} so that for every x G 

b{x, z) = ip{x^)ip{\z\) a.e. z€M.^, 

where (p : —)> M is a non-negative measurable function, and if : M+ —>• M is a measurable 

function such that 

if{r) 

, „ is non-increasing in r > 0. (1-7) 
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Theorem 1.5. Suppose Assumption 1 and Assumption 2 hold. Let H G : Xi> r(x*)| he 
a special Lipschitz domain, where F : —>■ M is a Lipschitz function with Lipschitz constant Aq 

(that is, |r(x*) — r(y*)| < Ao|x* — y*| for every G Then there are positive constants 

Ri = Ri{d,a, /3, Xq, Ml, M 2 ) and C 3 = C 3 {d,a, (d, Xq, Mi, M 2 ) > 1 such that for every r G (0, i2i], 
there is a constant rji = rii{d, a, (3, Xq, Mi, M 2 , r) G (0, r/2) so that for every zq G dD and every non¬ 
negative function f that is harmonic in DPiB{zQ, r) with respect to and vanishes in D‘^nB[zo, r), 

< |v/(x)| < forxeDnBizo,Vi). (1.8) 

Od{x) od[x) 

Obviously Assumption 2 is implied by 

Assumption 3. There exists a measurable function if : M+ —)■ M satisfying (HZD such that for 
every a; G 

b{x, z) ='ipdz]) a.e. 2 G 

Definition 1.6. An open set D C is said to be Lipschitz if for every zq G dD, there is a Lipschitz 
function F^g : —>• M, an orthonormal coordinate system CSzg and a constant Rzg > 0 such 

that if y = (yi, ■■ ■ , yd-i,yd) in CSzq coordinates, then 

DnB{zo,Rzo) = {y ■yd> Tz^iyi,--- ,yd-i)} n B{zo,Rzo)- 

If there exist positive constants Rq and Aq so that Rzq can be taken to be Rq for all zq G dD 
and the Lipschitz constants of F^g are not greater than Aq, we call D a Lipschitz open set with 
characteristics {Xo,Ro). 

Clearly, if D is a Lipschitz open set with characteristics {Xo,Ro), then it is K-fat for some 
K = fi;(Ao, Rq) S (0,1). The following theorem follows directly from Theorem 11.51 

Theorem 1.7. Let D be a Lipschitz open set in with characteristics {Xo,Ro). Under Assump¬ 
tions 1 and 3, there are positive constants R 2 = R 2 {d, a, (3, Xq, Rq, Mi, M 2 ) and 
C 4 = Ci{d,a, (3, Xq, Rq, Ml, M 2 ) > 1 such that for every r G (0, i? 2 ]; there is a constant 772 = 
r} 2 {d,a, (3, Xq, Rq, Ml, M 2 ,r) G (0, r/2) so that for every zq G dD and every non-negative function 
f that is harmonic in D n B{zQ,r) with respect to and vanishes in D'^ n B{zo,r), 

< \Vf{x)\<C^P^ forxGDnB{zo,m)- 

Onyx) od{x) 

Results in Theorem ll.41 Theorem 11.51 and Theorem 11.71 can be called uniform gradient estimates 
because the constants Ck, 2 < fc < 4, and rji, 1 < i < 2, are independent of eo of (|1.6p . 

The rest of the paper is organized as follows. Preliminary results on Green functions and Poisson 
kernels are presented in Section [2l The proof of the uniform boundary Harnack principle is given 
in Section [3l Section S] is devoted to the proof of Theorem 11.41 while the proof of Theorem 11.51 
is given in Section [5l In this paper, we use capital letters Ci,C 2 ,--- to denote constants in the 
statements of results. The lower case constants ci,C 2 , • • • , will denote the generic constants used 
in proofs, whose exact values are not important, and can change from one appearance to another. 
We use Ck to denote the unit vector along the positive direct of x^-axis. 
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2 Preliminaries 


Recall the Levy system [J^{x,y)dy,t) from (jl.ljl . which describes the jumps of X^\ for any non¬ 
negative measurable function / on M_|_ x x with f{s,y,y) = 0 for all y G x € and 
stopping time T (with respect to the filtration of X^), 


E. 


J^f(s,Xl,X,^) 


s<T 


= E. 


/ f{s,Xl,y),f{Xl,y)dyds 


( 2 . 1 ) 


Suppose II is a Greenian open set of X^. Let G\){x,y) denote the Green function of D, that 


IS, 


f[y)G\){x,y)dy = E^ 


ID 


/ f{X^s)ds 

Jo 


for every bounded measnrable function f on D and x G D. It follows from (12.ip that for every 
bounded open set D in every / > 0, and x G D, 


E, 




/ fiz) 

Jd<^ 


'D 


G'h{x,y)J^{y,z)dy dz. 


Defi n e 




y{x,z):=f GD{x,y)J^{y,z)dy ior {x, z) G D x D^. 
Jd 

We call K\){x, z) the Poisson kernel of X^ on D. Then (j2.2p can be written as 


( 2 . 2 ) 

(2.3) 


E. 


/(X(:j : + X 


I D<^ 


f{z)K^j^{x,z)dz. 


For any A > 0, define 


(2.4) 


h\{x,z) \= \d "6(A ^x,A ^z) forXjZGM'^. 

It is not hard to show that 

jG{x,y) = A-('^+")j^(A-^x, A-^y) for x,y G M'^. 

and 

> 0} has the same distribution as {X^-^t > 0}. 

So for any A > 0, we have the following scaling properties; 

G’hix, y) = X‘^~°‘GI)j{Xx, Xy) for x,y G D, 

K\){x, z) = X'^kGj^{Xx, Xz) ioT X G D, z G D'^. 

If u is harmonic in D with respect to X^, then for any A > 0, v{x) := u{x/X) is harmonic in XD 
with respect to X^^. 

When b{x, z) = 0, X^ is simply an isotropic symmetric a-stable process on which we will 
denote as X. We will also write J for J^. It is known that if d > a, the process X is transient and 
its Green function is given by 

^id/2) _Aa-d r _ ^ md 


(2.5) 

( 2 . 6 ) 


G{x,y) = 


2«7r'^/2r(a/2)2 


\x — y|“ for x,y G 


(2.7) 
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It is shown in Blumenthal et al. [T] that the Green function of X in a ball B{0,r) is given by 


GB( 0 ,r){x,y) = 2aJi/2Y{l/2y ^ ^ |x - y|“ foT x,y € B{0,r), (2.8) 

where z = (r^ — |xp)(r^ — \y\‘^)\x — y\~‘^ and r > 0. The above formula yields the following two-sided 
estimates (see, for example, my- Suppose B is an arbitrary ball in with radius r > 0. Then 
there is a universal constant ci = ci(d, a) > 1 so that for every x,y € B, 


GB{x,y) X 


x — 



6b{x) 

1 ^ - y\J 



^Bjy) 

- y\ 


a/2 


(2.9) 


Since for a, 6 > 0, o A 6 x and 1 A | in view of (11.31) we can rewrite (|2.9I) as 


GB{x,y) X |x - ?/|" 


5B{xr/^5B{yr/^ 

rB{x,y)°‘ 


( 2 . 10 ) 


It follows immediately from (|2.9I) that there is a positive constant C 2 = C 2 (d, a) > 1 so that for 
B = B(xo,r), 

Cg < 'RxTb < C 2 r" for x G B{xo,r/2). 


Riesz (see m) derived the following explicit formula for the Poisson kernel z) of X on 

B{0, r). 


KB{Q,r){x,z) = 


r(d/2) sin(7ra/2) (r^ — |xp)"/^ 


TT' 


■d/2+1 _ ^2^a/2 l^. _ 


for |x| < r and \z\ > r, 


( 2 . 11 ) 


We also know from [21 Lemma 6] that, for rotationally symmetric a-stable process X, 


^x{Xro G dD) = 0 and so Px(^ro / ^r^-) = 1 

for every x € D if D is a, domain that satisfies uniform exterior cone condition. (Although the result 
of [H Lemma 6] is stated only for bounded domain satisfying uniform exterior cone condition, its 
proof works for any domain satisfying uniform exterior cone condition.) 


3 Boundary Harnack principle 

Recall that we write X and J for X^ and J^, respectively. First we record the following gradient 
estimate on the Green function Gb oi symmetric a-stable process X from [3]. 

Lemma 3.1. ([H Gorollary 3.3]) Let D he a Greenian domain in of X. Then 

\VGD{x,y)\<d- — forx,y e D, x^y. (3.1) 

\x -y\ Ado(x) 
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For X ^ y in D, define 


hD{x,y) := < 


lx — y|" ^ (1 A 


— y\l^ ( 1 A 


^piy) 

\x-y\ 


- y\ 

The following two results are established in m- 


\^-y\) 

«/2 / I |\/3-a/2 


if a > 2/3, 
if a = 2/3, 


(3.2) 


1 1 A ■ (ivS) ' if a <2/3, 


Lemma 3.2. ([TSl Theorems 4.11 and Lemma 4.13]) Suppose b is a bounded function satisfying 
(ILSI), and that for every x G J^{x,y) > 0 a.e. y G M'^. There exist positive constants ri = 
ri(d, a,/3, Ml) G (0,1] and = C^{d,a, such that for any xq G and any ball B = 

B{xo,r) with radius r G (0,ri], we have for x,y ^ B, 

^GB{x,y) < GYx,y) < ^Gb(x,?/) and |5^G^(x,?/)| < G^hsix^y) (3.3) 

Moreover, Pa; {Yb ^ ~ every x G 33. In this case, for every non-negative measurable 

function f, 

^xf{xY=[ f{z)KB{x,z)dz forxeB. 

Jb^ 

Lemma 3.3. f [T3l Lemma 4.2]) Let D be a bounded open set in M'^. There exists a constant 
Ce = C' 6 (d, a,/3, diam(Zl), Ml) > 0 such that for any bounded function b satisfying (11.51) . and that 
for every x G J^{x, y) > 0 for a.e. y € we have 


G^Yx,y)<Ge\x-yr-<^ forx,yeD. 


Note that the constant C 7 below is independent of eq S [0,1] appeared in (|1.6I1 . 

Theorem 3.4 (Uniform Harnack inequality). Let ri G (0,1] be the constant in Lemma \3.‘A Under 
Assumption 1, there exists a constant Gy = C 7 (d, a,/3, Mi, M 2 ) > 1 such that for every xq G 
r G (0, ri], and every non-negative function u which is regular harmonic in B{xQ,r), we have 


sup u{y) < Gj inf u{y). 

y£B(xo,r/2) y&B{xo,r/2) 


Proof. Let m*(x) := Ea 


u(X^\ ) 

' ^B{xq , r ) ' 


Then u* is regular harmonic in B{xo,r) with respect to 
land Assumption 1, for every xq G and 
r G (0, ri], the Poisson kernel Kb{xo B{xo,r) of is comparable to that of Thus 

for every x G B{xo,r/2), u{x) is comparable to u*{x). Theorem 13.41 then follows from the uniform 
Harnack inequality for mixed stable processes; see [3 (3.40)]. □ 


the mixed stable processes X^°. In view of Lemma 


Lemma 3.5 (Harnack inequality). Under Assumption 1, there exists a constant Gs = G^{d,a, jl, 
Ml, M 2 ) > 0 such that the following statement is true: If xi,X 2 G r G (0, ri] and k are 
such that |xi — X 2 I < 2 ^r, then for every non-negative function u which is harmonic with respect to 
X^ in B{xi,r) U B (x 2 ,r), we have 

^_l2-A;(d+a)^(x2) < y^^xi) < C82^('^+“)u(x2). 


(3.4) 











Proof. Without loss of generality, we may assume |xi — X 2 I > 'r/4. Note that for every x € 
B{x 2 ,r/ 8 ) C B(xi,r/ 8 )‘^, we have \x — xi\ < Thus by Lemma [32] and Assumption 1, we 

have 

^ [ GB{xi,r/8){xi,y)J{y,x)dy 

JBixi,r/8) 

~ 2 M 2 x/8) (xijX) 

> = C2r-'^2-'=('^+“). (3.5) 

Recall that by Theorem 13.41 we have u{x) > c^u{x 2 ) for every x G B{x 2 ,r/ 8 ). Thus by (13.51) . 


n(xi) > / u{x)KL ,^Jxi,x)dx 

Jb{x 2 x/ 8 ) 

> C 2 C 3 n(x 2 )r-'^ 2 -^('^+") [ dx 

J B{x 2 ,r/ 8 ) 

> C42-’^G+<^^u{x2), 


and (13.41) follows by symmetry. □ 

Proof of Theorem 11.31 Note that there are constants Rq = Ro{d,a, (3, Mi) G (0,ri) and c = 
c{d,a, P, Ml) > 1 so that 

- - rr— < J^{x,y) < - -ri— 

\x-y\d-+a- ^ ^^ - \x - y\d.+» 

for all \x — y\ < Rq and b(x,z) satisfying (|1.5p . Thus using (|3.3p . we can get uniform estimates on 
the Poisson kernel 

■J B{xQ,r) 

of any ball B{xQ,r) with respect to with r G (0,i2o/3), x G B{xo,r) and r < |z — xo| < 2r. 
Specifically, for r < \z — xq\ < 2r, K^ixo uniformly comparable to KB(xo,r)i^^ Using 

the explicit formula ()2.11l) for the Poisson kernel Kb^xqx)^ (13.3p . Theorem 13.41 and (II.6 p . we can 
adapt the arguments in jS] Theorem 2.6] to get our uniform boundary Harnack principle 11.31 (cf. 
the proof of [T] Theorem 3.9]). Since the proof is almost identical to those in [S] Section 3], we omit 
the details here. □ 


Lemma 3.6. Suppose Assumption 1 holds and D is a Lipschitz open set with characteristics 
(Ao,Ro)- Let ri G (0,1] be the constant in Lemma [g.gl There is a positive constant Cg = 
CQ{d, a, j3, Xq, Rq, Ml, M 2 ) > 1 such that for every zq G dD, r G (0,ri/2), and every non-negative 
harmonic function u that is regular harmonic in D Pi B{zo, 2r) with respect to and vanishes in 
n B{zo,2r), 


E, 




b 

TDCBi^ 




<Cg2-‘^^u{x), x^DPBk 


for Bk '■= B{zq, 2 ^r) and k>\. 


(3.6) 
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Proof. Without loss of generality, we may assume zq = 0. By the uniform inner cone property of 
a Lipschitz open set, one can hnd a point zq G D n B{0,r) and k = k{Xo,Ro) G (0,1) such that 
Bjf := B{zk, n2~^r) C B^Ci D for every z^ ■= 2 ~^zq and fc > 0. Define 


Ufc(x):=E, 




Since uq = u, (IMI) is clearly true for k = 0. Henceforth we assume k > 1. Note that rtfc > 0 is 
regular harmonic with respect to X^ in D D B^, and Uk{x) < Uk-i{x) for all x G Dehne 


4(x):=E, u{Xl^-.Xl^^ GB, 


Clearly by dehnition Uk{zk) ^ Ik{zk)- For any ^ > 1; by Lemma 13.21 and (|2.5I) . we have 

Ks^Zk^y) = [ G^Bkizk,z)J^{z,y)dz 

Jbu 


< - 


2 JB, 


GBk{h,z)j’'{z,y)dz 

G^-ik-i)B^ {2~^''~^'^zi,z)J^{z, y)dz 


2 72-('=-i)Si 

^ JSi 


^2 (fc i)q f GBiih,w)j’’{2 ^'>w,y)dw. 
2 Jbi 


Note that for any y G Bq and w G Bi, 

\y - w\ 


< 


\y\ + kl 


< 3. 


\y — 2~^^~^^w\ \y\ — 2 ^+^|u)| 

Thus by (11.61) we have 

k(2-(^-bu;,y) < M2J"°(|y-2"^^"^k|) < 3^+^M2r°{\y-w\) < 3^+^Mlj\w,y). 
It follows from (j3.7l) . (13.8|) and Lemma [22] that for any y G Bq, 

od+CK+l r 

Kski^k.y) < —5—M|2“(''"b« / GBA2i,w)J^{w,y)dw 

2 JBi 


< 3'='+“+iM|2-('=-b“ [ G’’B,ih,w)j\w,y)dw 

Jbi 


= ci2-^“iL^^(5i,y). 


Now we have for A; > 1 


4(4) = / u{y)KBi^{zk,y)dy 

Jb- 

< ci2“^“ [ u{y)K^Bii2i,y)dy 

Jbs 


(3.7) 


(3.8) 
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(3.9) 


Next we compare Ii{zi) with u{zi). Using Lemma [321 (|1-6I) and (12.5p . we have 


= f 

J u 


\z—zi\<Kr/2 


G^~^{zuz)j\z,y)dz 


ki 

ki 


> - / Gg (zi, z)J^°{\y — z\)dz 


Gi^BiiO,z- zi),r°{\y - z\)dz 

\z—zi\<K,r/2 

[ G^Bi{0,w)J^°{\y - zi-w\)dw 

J\w\<Kr/2 

f G^Bi (0, Kz)J^°{\y -h- Kz\)dz 
J\z\<r/2 


1 


/ GBiiO,z),P°{\y - zi - Kz\)dz. 


2M2 JBi 
Again using Lemma 13.21 and (|1.6|) , we have 


= / G%i{zi,z)r{z,y)dz 

JBi 

< ^M 2 [ GBiih,z)J^°{\y - z\)dz 

z Jbi 


= Im 2 


= 1^2 


+ 


'LI<Ll|/2 ./|zi|/2<|z|<r/2 

GBiih,z)J^°{\y - z\)dz 


GBA^iiz)J^°i\y - ^\)dz 


'kl<Lil/2 


+ 2 “ 


' |zi|/4<|?ii+ii/2|<r/4 


Gbi{zi,2w + zi)J^°{\y — zi - 2w\)dw 


Note that for any y G and \z\ < |2i|/2, |z — zil > l^i] — 1^1 > 1^1 and \y — zi — Kz\/\y 
{\y\ + l^il + I^k\) / (|y| - \z\) < 4. Thus 


GBiizi,z) X | 2 - 2 :ir lA 


a-d I 


Iz-zil^ 


< u 




and 


r°{\y - ^1) < r^{\\y -zi- KZ\) < 4''+“J^“(|y - - KZ\ 


It follows then that for any y G 


0 ’ 


[ Gsi(2i, 2;)J^°(|?/- z|)(iz < C2 / GBiiO,z)J^°{\y - zi - Kz\)dz. 

■^L|<|ii |/2 J\z\<\z ^\/2 


LI<|A|/2 


(3.10) 


(3.11) 

- 2I < 


(3.12) 
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Note that for y £ Bq and |5i|/4 < \w + zi/2\ < r/4, 6bi{‘2‘W + zi) = r/2— \2w + zi \ < 2(r/2 — |t(;|) = 
25 biM, and \y-zi-Kw\/\y-zi-2w\ < (|y| +k\w + zi/2 \ + (1 - k/2)\zi\) / (|y| - |zi +2w\) < 2. 
Thns 


GBiih, 2 w + h) - 


< 


\2w 


a—d 


1 A 


m; 


\a—d 


1 A 


Sbi ( 2 w + 

| 2 t(;|“ 


U! ' 


(0,u;), 


and 


r°{\y -zi- 2w\) < r°{\y - zi - ku;|/2) < 2'^+"J^°(|y - - kw\). 

Thns for any y £ 


' \ z\\/ 4:<\w+z\ /2\<r / A 


Gbi{zi,2w + zi).P°{\y — zi - 2w\)dw 


< C 3 [ GBiiO,w)J'^°{\y - zi - Kw\)dw 

J\zi\/A<\w+zi/2\<r/A 

< C 3 / GBiiO,w)J'^°{\y - zi - Kw\)dw. 

Jbi 

Using ()3.12l) and (I3.13p . we can continne the estimates in (|3.11l) to get that for any y £ Bq 


K' 


Bi{h,y)<C4 [ - Zi - Kz|)d, 

JBi 


z 


(3.13) 


(3.14) 


Combining ()3.10p and ()3.14p . we get 


K%(h,y)<C5K'^K^~(zi,y), foryG.B| 


Bi' 


It follows that 


h{h) = [ u{y)K^B^{zi,y)dy <C 3 K [ u{y)K^~ {zi,y)dy 

Jbs Jb^ "1 

< csK"" / u{y)K^~^{zi,y)dy = C 5 K~‘^u{zi). 

J 

Conseqnently by (|3.9I) and (I3.15P we have for all A: > 1, 


(3.15) 


Ukih) < Ik{zk) < C1C5K "2 ^“n(2i). 


(3.16) 


By the monotonicity of in k, Theorem 11.31 (|3.16l) and Lemma 13.51 we conclude that for any 
X £ D f] Bk and k > 1 


Ukjx) ^ nfc-i(a:) nfc-i(4-i) 
u{x) ~ u(x) ~ u(zk-i) 


< 


C6CiC5K-“2-(^-l)" 


^(^l) 

u(zk-i) 


< C 72 -^“. 


The proof is now complete. 


□ 
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For a Lipschitz open set D with characteristics (Aq, Rq), let k = «:(Ao, Rq) G (0,1) so that D is 
K-fat. For zq £ dD and r £ (0,1], we use Ar(zo) to denote a point in D such that B{Ar{zo), nr) C 
D n B{zo,r). 

The following lemma follows from Theorem 11.31 and Lemma 13.61 (in place of [21 Lemma 13 and 
Lemma 14]) in the same way as for the case of symmetric a-stable process in [21 Lemma 16]. 
We omit its proof here. See also HI Lemma 3.4] for detailed computation of a similar result for 
symmetric Levy processes. 

Lemma 3.7. Suppose Assumption 1 holds and D is a Lipschitz open set with characteristics 
(Ao,i?o)- Let ri £ (0,1] he the constant in Lemma \3.‘A There exist positive constants 71 = 
'fi{d,a, 13, Xq, Rq, Ml, M 2 ) and Cio = Cio{d,a, j3, Xq, Rq, Mi, M 2 ) such that for every zq £ dD, 
r £ (0,ri/2) and all non-negative functions u, v that are regular harmonic in D D B(zQ,2r) and 
vanish in n B{zo,2r) with u{Ar{zo)) = v{Ar{zo)) > 0, we have 

(i) h{zo) := \im.D^x-,zo u{x)/v{x) exists; 

(ii) 1^ -^(^o)| < Cio forx £ Dr\B{zo,r). 


4 Gradient upper bound estimates 

We now study gradient estimates for non-negative harmonic functions of in open sets. 

Lemma 4.1. Suppose b is a bounded function satisfying (II3I), and that for every x £ L^(x, y) > 

0 a.e. y £ M'^. Let ri £ (0,1] be the constant in Lemma and B = B{xo,r) with r £ (0,ri]. 
Then for every x £ B, z £ B^ and I <i < d, 


dxi [ GB{x,y)j\y,z)dy = [ dxiGB{x,y)j’’{y, z)dy, 
Jb Jb 


(4.1) 


dx. 


and 


GB{x,y)SlG^B{y^w)dy ) J\w,z)dw = 


B \Jb 


B \JB 


dxiGBix, y)SyG%{y, w)dyj J\w, z)dw, 

(4.2) 


dxiK%{x,z) = [ dxiGBix,y)J^iy,z)dy + 
Jb 


dxiGBix, y)SyGB{y, w)dy ) J%w, z)dw. (4.3) 


IB \JB 

Proof. Without loss of generality we assume i = d. Fix x £ B and 2 £ B'^. We have 

,A{d,-/3) 


I Tbi M / A[d,—a) |i^iicxD- 

Thus (|4.ip follows directly from Lemma 5.2]. 

Let Qziy) ■= jQ<SyG^^{y,w)J^{w,z)dw for y £ B. We have 


< - 1 - 00 . 


dxd / GB{x,y)gz{y)dy =\\m. 

Jb ^^oJb 


Gb{x + Xed,y) - GB{x,y) 
X 


9z{y)dy. 


(4.4) 
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To prove (14.21) , we only need to show that the integrand in the right hand side of (I4.4p is uniformly 
integrable on S in A G ( 0 , (5b(x)/ 2). Note that we have 

GB{x,y) = G{x,y) -Ey [G{x,Xrg)] =: G{x,y) - H{x,y). 


Thus 

\Gb{x + Xed,y) -GBix,y)\ ^ \G{x + Xed,y) - G{x,y)\ \H{x + Xed,y) - H{x,y)\ 
X - A A 

= : I+ 11. 


Obviously by ( 12 .7|) we have 

/ < Cl (|x + Acrf - + |x - 2/1“-'"-^) for y G B, 

for some positive constant ci = ci {d, a) . Since H{x,y) =Ey [G{x, )]) by the mean-value theorem, 

there is a point x\ in the line segment connecting x with x + Xed so that 

11 = d^^H{xx,y) = Ey[d^^G{xx,Xrg)] < C25 b{x)’^~^~'^■ 


Thus for some positive constant C 3 = C 3 {d,a,x), we have 

\GBix +Xed,y) - GB{x,y)\ ^ , x i., ixa-d-i , ,, 

-- < ci{\x + Xed-y\ A\x -y\) -hcs. (4.5) 

Let h{y) := hB{y,w)dw for y G B. Note that by Lemma [.S.2I and the boundedness of tc i->- 
J^{w, z) on B, 

\gz{y)\<CA j hBiy,w)j’’{w,z)dw < C5h{y). (4.6) 

Jb 

Thus by (14.5p and (14.6p the integrand in the right hand side of (14.41) is uniformly integrable on B 
in A G (0,(5 b(x)/ 2) if the following three conditions are true: 


(i) !b Ky)dy < + 00 ; 

(ii) su-p.^^b{x,Sb{x)/ 2 ) JB^^y)\y - dy < -Foo; 

(hi) lim^io SUp^GB(x,5sW/2) J{yeB:\y-w\<e} d{y)\y “ ^ 

If a > 2/3, then for any y ^ B, 


Kv) < [ \y- w\‘^-^-'^dw < [ \u\'^-d-<^du < + 00 ; 

J wGB J\u\<2r 

that is, h{y) is bounded from above on B. Obviously (i)-(iii) hold for h. If a = 2/3, then 

5B{w)d 




I. 


\w 


w^B,\w—y\>e5B{y) 


|/3/2-d dBjwf f 5B{yfG 

5B{y)dl‘^ dB{y) j 


dw 
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< 

r\j 


+ 


/ 

Jwi 


'w€B,\w-y\<eSB(y) 

+ 1 - 


|t(; — yf '^dw 


(4.7) 


Using this upper bound, it is easy to check h satisfies (i) and (ii). As for (hi), note that 6b{w) > 
5b{x)/2 for every w G B{x,5b{x)/2). Consider an arbitrary e G {0,6b{x)/4:). Then B{w,s) C B 
and 5s(y) > (5 b(x)/ 4 for every y G B{w,e). We have by ()4.7I) . 


'\y—w\<e 


Hy)\y - w 


\a—l—d 


dy 


< 


< 

rs_/ 


{6B{y)-^/^ + l)\w-yr-^-<^dy 

'\y-w\<e 

[ {6B{x)-f^/^ + l)\w-yr-^-%. 

J\y—w\<e 


Thus condition (iii) is implied by the fact that 


lim sup [ (SBix)-^/^+ l)\w-y\'^-^-'^dy = 0. 

w£B{x^Sb{x)/ 2) J\y—w\<£ 


When a < 2/3, similar to (|4.7p we have 


h{y) X / \w-yf 
J w^B 



5b{wY \ 
\y-w\^ ) 


/ \y-wf 

dBiyf-^/^ ) 


dw < dB{y)^ + 1 . 


By a similar calculations as in the case a = 2/3, we can show that (i)-(iii) hold for h. This completes 
the proof. □ 


Lemma 4.2. Under Assumption 1, there exists a constant Cn = Cu{d, a, (3, Mi, M 2 ) > 0 such 
that for every B = B{xq, 1) and 1 < i < d, 


L 


\dxiGB{x,y)\J^{y,z)dy <Cii GB{xo,y)J^{y, z)dy for x € B{xo, l/A) and z G B^. (4.8) 


IB 


Proof. Without loss of generality, we assume xq = 0 and i = d. For every \x\ < 1/4 and \y\ < 1, 
we have \x — y\ /\ 5b{x) x |x — y\. Thus by (13.ip . 


/ \dx^GB{x,y)\J^{y,z)dy 
J B 


< 


GB{x,y) hr 

J {y,z)dy 


'\y\<l \x-y\ AdBix) 
^B{x,y 
^\y\<l 1^ “ y\ 


GB{x,y) hr ^ , 
J {y,z)dy 


+ 


\Jl/2<\y\<l J\y\<l/2j \x - y\ 

= : I{x, z) + II{x, z). 


Ay.z)iy 


For 1/2 < |y| < 1 and \x\ < 1/4, we have |x — y| x \y\ x 1 and 6b{x) x 1 . Thus 


GB{x,y) X |y|“ 



<^s(2/)"^^\ 


GBif3,y), 
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and consequently 


I{x,z) X [ 

hi 


(4.9) 


G'B(0,y)J^(y,z)dy. 

For every \y\ < 1/2, |x| < 1/4 and \z\ > 1, we have \z — y\ > 1/2, \z — y\ x \z — y + x\ and 
6B{y) = 1 - |y| X 1 - |y - x| X 1. Thus by ([IT]) 

5B{x)h‘^5B{y)h^ 


II{x,z) 


ci(M2) 


<\y\<l/2 


|a: - y\ 


/ \x-y\ 

J\y\<l/2 


— d+OL—1 


—d+a—1 


1 A 


1 A 


lx - ?/|" 

(5b(x)"/^(1 — \y — 


\w 


\x-y\ 

-d+a-l ^ ( 1 - 


h°{\z-y\)dy 

h°{\z-y + x\)dy 


'\w+x\<l/2 




< 


C2(M2) 


L 


\w 


—d-\-OL—l 


1 A 


(1 — \w\)°'h 


|ui|<3/4 


U) ' 


h°{\z — w\)dw 


J^°{\z — w\)dw 


\w 


—d+a 


-1 Tb 


J\w\<^/A 

=■■ gi{z)- 


U) 


J°{w, z)dw 


For every z > 1, let 


92{z) := f 

J k 


1 A 


(l_|y;|)«/ 2 ' 


|ui|<3/4 




h{w, z)dw. 


Obviously 


( 72 ( 2 ;) X f GB{0,w)h{w, z)dw < I GB{0,w)h{w, z)dw. (4-10) 

j\w\<3/4: Jb 


' |ki|<3/4 

Note that J^°(|x — y\) is non-increasing in |x — y\. Thus by (ll.Op 


sup , . y I l\ /O \ 

Id>i 92iz) |,|>1 M^-V-od^l + 3/4) (l A ) dw 

where M = M{d, a, /3, M 2 ) > 0. Thus by (I4.10p and ()4.1ip we prove that 


giiz) M 2 hh\z\ - 3/4) |4n|“ ^ (l A dw 


< sup 


< M < -|- 00 , (4:.11) 


c3(Af2) r 

II{x,z) < / GB{0,w)J°{z,w)dw for |x| < 1/4, |z| > 1. (4-12) 

Jb 

Therefore (14.81) follows from ()4.9p and ()4.12p . □ 

Recall the definition of r£){x,y) and h£){x,y) from (11.21) and (18.21) . respectively. 

Lemma 4.3. For B = i?(0,1), there exists a constant G 12 = Gi 2 {d,a, fd) > 0 such that for every 
1 <i <d, |x| < 1/4 and |r(;| < 1, 

[ \dxiGB{x,y)\hB{y,w)dy < Gulx - _ 

Jb 
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Proof. Without loss of generality, we assume i = d. For any |x| < 1/4 and \y\ < 1, we have 
5b{x) X 1 and \x — y\ /\ 5b{x) x |x — y\. Thus by (j3.ip . 


'B 


\dxdGB{x,y)\hBiy,w)dy 


< 


GBix,y) 


'|y|<i \x-y\/\ Sb{x) 

f GBix,y) 

'\y\<i \x-y\ 


hBiy,w)dy 


hB{y,w)dy. 


(4.13) 


We note that for |x| < 1/4 and |y| < 1, 


GB{x,y) X |x - ?/|" I 1 A 


_y\a/2 


1 A 


Ssiy) 


a 12 


_y|a/2 


\X 


\x - y\ 


rB{x,yY 


(4.14) 


and rB{x,y) > Sb{x) > 3/4. Now we calculate the integral in (j4.13p using (j4.14p and the explicit 
formula of hB{y,w). If a > 2(3, we have 


(IM) 


(5b (re) 


a/2 


^B{y) 


a/2 


J\y\<i ^+1^ rB{x,y)<^rB{y,w)°‘G 

< <5b(^u)“/2 ' 


dy 


< 


\x — w 


\x_y^-d+a-l\y_^^-d+a-ydy 

-d+{a-l)h{a-0) (y;)"/2 


< |x-u;|-'='+("-^)''("-^)<5b('«;)^. 

If q; = 2/3, we have by (|2.10l) and (13.21) 

(imi) 

rB{y,w)P V 


(4.15) 


/ \x-y\ 

J\y\<^ 


2/3 


i d dB{yf 


rB{x,y)^l^ 

5b{wY 


^B{y)^ 


+ 


\y\<i,\y-w\<(^&B{y) \x - ‘^f^\y - w\<^ d rB{x,y)^drB{y,w)d 

6b{w)I^ 


^B{y) 

dy 


dy 


\y - w\f^G§^(^y)y/‘i I 6B{y)^G |y - 'w| \ 


l\y\<l,\y-w\>e5Biy)\x 2/3|y;_y|d /3 y;)/3 |y _ y;|/3/2 ^^(y) 


< 

rs_/ 


5b{wY / ^ ‘^\y-wf '^dy 

J \y\<^,\y-M<(^^B{y) 


+6b{w)^ / \x — y^^ ^ ‘^\y — ’^dy 

J\y\<^/y-w\>^^B{y) 

< 5B{wf\x - = |x - 


(4.16) 


If a < 2/3, we have 


(iraii 


'iyi<i 


\x - y\~ 


< 


_rf+„_i I _ ..A-d+a -0 5b{wTG / \y-w\P 

rBiy,w)'^G ^ 5B{yY~°'G 

I'w^ “ y\^~°'G5B{y)^~^ 


I \a\y~'^ 

^B\x,yr 

5b{wTG 


dy 


'\y\<l,\y-w\>5B{y) \x - y\^ °‘+^\w-y\^ "+^ rB{y,w)<^G 


dy 
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+ 


6b{w) 


a 12 


^bW) 


a 12 


\y\<i,\y-w\<SB(y) \^ - - y\'^ "+/^ rB(a;, ?/)"rB(?/, u;)“/2 


dy 


< 

rs_/ 


SB{wr^^ [ 

Jk 


\y\<l,\y-w\>5B(y) 


_y\-d+a-l\y_w\-d+»-Pdy 


J\y\<^^y-^\<dBiy) 

Sb{wT^‘^ [ \x-y\-^+^-^\y-w\-^+^-^dy 

J\y\<i 


< 6b{w)^^‘^\x - y;|-<i+(«-l)A(a-/3) f (1^ _ y|a-l-(«-l)A(«-/3) ^ _ ^|a-/3-(a-l)A(a-/3) 

■>'l!/|<l 

< |x-(4.17) 

Lemma 03] follows from (I4.16p . (I4.17P and (I4.15p . □ 

Lemma 4.4. Under Assumption 1, there exists a constant C 13 = Ci 3 {d, a, j3, Mi, M 2 ) > 0 such 
that for B = B{xq, 1), 1 < i < d, x ^ B{xq, 1/4) and z G B'^, 


iB 


IB 


\dcciGB{x,y)S^G^Biy,w)\dy 


J^{w, z)dw < Gi 3 / Gb{xo,w)J^{w, z)dw. 


IB 


(4.18) 


Proof. Without loss of generality, we assume that xq = 0 and i = d. Let ri G (0,1] be the constant 
in Lemma 13.21 By Lemma 13.21 and the scaling property, we have for y,w (z B, 


\StGUy,w)\=rf\Sl’'^ 


G'r?s(ny,n«^)| < cirf\h ri B{riy,riw)\ = cir“ ^hB{y,w) < cihB{y,w). 


Here ci = ci{d,a,l3, Mi) > 0. Hence to prove (14.1811 . it suffices to prove that for x G H(0,1/4) and 
z G B^, 


B 


B 


\dxaGB{x,y)\hB{y,w)dy 


J^{w,z)dw<C 2 / Gb{xo,w)J^{w, z)dw (4-19) 
Jb 


for some C 2 = C 2 {d, a, f5, Mi, M 2 ) > 0. By Lemma 031 we have 


< 

rsj 


B UB 


\dxdGBix,y)\hBiy,w)dy 


J^{w, z)dw 


J |to|<l 


+ 


dB{w)^^^\w - J'’(u;,z)du; 


4|to|<l/2 4l/2<|to|<l 

= : I{x, z) + II{x, z). 

Fix |x| < 1/4 and \z\ > 1. For 1/2 < |tc| < 1, we have |tc — x| x |u;| x 1, and consequently 
Gb{0,w) X (5s(r(;)“/^. Thus 

//(x,z)x f 6 b{w)°'^‘^J^{ w,z)dw ^ f GB{^iw)J^{w,z)dw. (4.20) 

Jl/2<\w\<l Jl/2<\w\<l 
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For any |u;| < 1/2, we have 1/2 < 6b{w) < I, \z — w\ > 1/2 and \z — w + x\ x \z — w\. Hence by 

m) 


I{x,z) 


Ma 


_ |-(i+(o-l)A(o—/3) 


\W — X 


'\w\<l/2 


J^°{\z — w\)dw 


[ l-u; - J"°(|z -w + x\)dw 


J\v+x\<l/2 


< 


■lhl<3/4 


(4.21) 


92 


(4.23) 


We first consider the case \z\ > 2. Let g 2 iz) := ii„|< 3/4 '^J^°{\z — v\)dv. Note that for any 

< 3/4, we have Gb{0,v) >; Thus 

(z) X f GBiO,v)J'^°{\z — v\)dv < M 2 f GB{lliv)J^{v,z)dv. (4.22) 

4hl<3/4 4|^|<3/4 

In addition since J''°(|?/|) is non-increasing in |?/|, we have 

Si(A ^ - 3/4) 4,^3,^ |„|-4+(„-l)A(„-«rf„ 

sup -^ < sup ----- - 3 --- < M < -1-00, 

\z\>2 92 { z ) |^|>2 + 3/4)/|^l<3/4|r;|" 

where M = M{d, a, /3) >0. Therefore by (14.211) . (14.221) and (|4.23l) we have 

csfMa) r 

I{x,z) < / Gb{11,w)J°{w, z)dw for Ixj < 1/4 and \z\ > 2. 

Jb 

On the other hand if 1 < \z\ < 2, we have 0 < (5 b(-s) < 1, and by (|2.1ip 

f GBiO,w)J^{w,z)dw > M2^ f GB{0,'w)J{w,z)dw = M2^KBiO,z) ^ M2^5Biz)~'^^‘^ > M2^- 
Jb Jb 

(4.25) 

Note that |z — rc| > 1/4 for any |rc| < 3/4. Thus 


(4.24) 


gi{z) < J"°(l/4) / < 1. 

4|«)|<3/4 


Thus by (I4.2ip . (I4.25P and p4.26p we have 

C4(M2) 


C4(M2j r 

I{x-,z) < / GB{0,w)J°{w,z)dw for |x| < 1/4 and 

Jb 


1 <\z\< 2. 


(4.26) 

(4.27) 

□ 


Now (I4.19P follows from (I4.27P and (I4.20p . 

Theorem 4.5. Let ri € (0, 1] be the constant in Lemma Id.21 Under Assumption 1, there exists 
a constant G 14 = < 714 ( 0 ?, a, /3, Mi, M 2 ) > 0 such that for every ball Br = B{xo,r) with radius 
r G (0, ri] and 1 < i < d, 


|5i,.iL5^(x,z)| < ^iL 5 ^(xo, 2 ;), for X € B{xo,r/A) and z £ Br 


(4.28) 
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Proof. Let \ := 1/r > 1/ri > 1 and define h\{x,z) = x, z). Observe that 1|6a||oo = 

r"“^||6||oo < ^Ml < Mi. By the scaling properties (j2.4jl and (j2.6jl . b\{x, z) satisfies Assumption 

1 and it suffices to show that for the ball B = B{xq, 1), 

IdxiK^^gix, z)\ < Ci 4 K^q{xo, z) for x € B{xo, 1/4) and 2 G B^. (4.29) 


We know from m Lemma 4.10] that 

G''B{x,y) = GB{x,y) + f GB{x,z)SGG^^{z,y)dz ioix,yeB. 

Jb 

Thus by (12.31) . for i = 1, • • • ,d, every x £ B, and 2 : G B'^, 


dxiK^B{x,z) = dxi j^GB{x,y)jG{y,z)dy+ dxi (^j^GBix,y)SGG^^{y,w)dy'^ jG{w,z)dw. 
Thus by Lemma l4.ll 

\9xiKGi^^^)\ ^ / \9xiGB{x,y)\jG{y,z)dy 

Jb 

+ Jb (^J^\dxiGBix,y)SGG’’jf{y,w)\dy'^ jG{w,z)dw. 

On the other hand, by (13.3p and (12.5p . we have 

^GB{x,y) < G^^{x,y) < ^GB(x,y) for x,y £ B. 

Thus 

K^B''ix,z)= f G^jf{x,y)jG{y^z)dy>^f GB{x,y)jG{y, z)dy for 2 ; G BL 
Now (|4.29p is implied by (13.3p and Lemmas I4.2M.41 This completes the proof of the theorem. □ 

Lemma 4.6. Suppose Assumption 1 holds and f is regular harmonic with respect to in B{x,r) 
for some x G and r £ (0,ri]. Then dx^fix) exists for every 1 < i < d and 

dxjix)= f _ J{z)dxiK'’B{x,r)ix,z)dz. (4.30) 

J B(x,r) 


Proof. Recall that is the unit vector along the positive Xj-axis. Choose e > 0 sufficiently small 
so that X + eei £ B(x,r/4). By the regular harmonicity of /, we have 

f(x + £ei) - fix) ^ r 

^ J B{x,r)‘^ 

Therefore (j4.30p follows from (|4.28p and the dominated convergence theorem. □ 


K{x,r)i^ + - ^B{x,r)(^^^) 


dz. 


Proof of Theorem 11.41 Let X £ D and 0 < r < (Sb(x) Ari)/2. Note that under our assumption, 
/ is regular harmonic in B(x,r) with respect to X^. By M.30I) and M.28I) . we have 


ldxj(x)l < 


' B{x,r) 


f{z)\dxiKB(^x,r) 


(x, z)\dz 
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< 


^ J B{x,r) 

Cm,, 

- fix) —^ ^ 

r ri A 6d{x 


[_ _ Ji^)KB{x,r)ix,z)dz 

J B(x.r) 

2Ci4 


f{x) as r t (ri A 5 d(x))/2. 


□ 


5 Gradient lower bound estimate 


For X = (xi, • • • , Xd) G we write x = (x, x^), where x = (xi, • • • , Xrf_i). In this section, we fix 
a Lipschitz function F : — )■ M with Lipschitz constant Aq so that |r(x) — F(y)| < Ao|x — y\ for 

all x,y € Put p{x) := Xd — F(x). Unless stated otherwise, D denotes the special Lipschitz 

open set defined by iA = {x G : p{x) > 0}. When x G iA, p{x) serves as the vertical distance 
from X G AA to dD, and it satisfies 

p{^)/\l 1 + -^0 — doix) < p{x) for X € D. (5.1) 

We define the “box” D~^{x,h,r) := {y G : 0 < p{y) < h, |x — y| < r}, and the “inverted box” 
D~ (x, h, r) := {y G : —h < p{y) <0, |x — y| < r}, where x G and h, r > 0. 


Lemma 5.1. Let ri G (0, 1] be the constant in Lemma [23 Suppose Assumption 1 holds, zq G dD 
and r G (0,ri/2]. Let £ D be such that p{Azq) = \Azq — zo\ = r/2. Then there exist positive 
constants C 15 = Ci 5 (d, a,/3, Aq, Mi, M 2 ) and 72 = 72 (d, a,/3, Aq, Mi, M 2 ) such that for every non¬ 
negative function u which is harmonic in D f] B{zq, 2r) and vanishes in D'^ D B{zo, 2r), we have 


u{x) ^ ( pjx) 

u{A,,)- ^^\p{A,,)J 


for X G iA n B{zo, r). 


Proof. Note that by Lemma l,'I.2l and Assumption 1, we have for every x G and y G B{x,r) 


K 


B{x,r) 


ix,y) 


M 2 


L 


B(x,r) 


GBix,r)ix,z)J^°iz,y)dz 


^B{x,r)(^^y)- 


(5.2) 


Lemma O follows from (j5.2|) . the uniform Harnack inequality (Theorem 13.4p and a standard 
argument of induction in the same way as for the case of symmetric a-stable process in [21 Lemma 
5] (see also [H Lemma 4.2]). We omit the details here. □ 

In the remaining of this section, we assume that Assumption 1 and Assumption 2 with i = d 
hold. In this case, the jumping kernel J^(x, y) of the process satisfies that for every x G 


J\x,y) 


Ad,-a) , p{x)'ijj{\y - x\) 

P) \y_^\d+0 


=■ j\x, \y 


x|) a.e. y G 


(5.3) 


We note that by condition (jl.7p of Assumption 2, j^(x, |z|) is non-increasing in \z\ for every x G 
g g (0, ri]. We define iA+ := lA+(zo, 4x^/1 -|- Aq, 2r) D IA n ^(zo, 2r), and 

gb,r{x) := Px ^ D~{zo, 00 , 2r)^ for x G M'^. 


Clearly, is regular harmonic in IA+ with respect to X^, gb,r{x) = 0 in IA {zo,oo,2r), and 
gb,r{x) = 1 in (IA+ U IA“( 2 ;o, 00 , 2r))''. 


21 












Lemma 5.2. The function gb^r{x) is non-decreasing in Xd- 

Proof. Note that gb,r{x) = 1 — Px ^ D~{zo, oo, 2r)^ for every x € Take x, y € such 

that X = y and yd < Xd- Consider the process {X^,Fx) starting from x (i.e. Px(^o = x) = 1). For 
every t > 0, define 

y/ := - {xd - yd)ed- 

Then {Y^,Fx) is a Markov process starting from y. Let denote the totality of tempered 

functions on For every / G if we define fd{z) := f{z — {xd — ydf^d) for z G then 


d"f{x - {xd - yd)ed) = C'fdix) = lirn f {fd{y) - fd{x))j''{x, \y - x\)dy for x G 

‘^^0J\y-x\>e 


Thus 


/(y/) -/(To') - f C^f{Y,^)ds = fd{X^)-fd{X^o)- fc^fd{X^s)ds 
Jo Jo 

is a Px-martingale. We know from m Theorem 1.3] that the solution of the martingale problem 
{C^jSfU.'^)) with initial value y is unique. Hence (Y^^,Px) has the same distribution as {X^^Fy). 
Consider the trajectory ui of X^ starting from x. If ui exits by going into +oo, 2r), then 

so does oj — {xd — yd)&d which is the trajectory of Y^ starting from y. Hence 


'Pa: (^T^+ e D {zQ,QO,2r)^ < Px e D ( 2 : 0 , 00 ,2r)^ = P^ G D {zo,co,2r)^ . 


This completes the proof. □ 

Lemma 5.3. Let ri G (0,1] be the constant in Lemm,a, \S. HI There are constants 

Cie = Ci6(d, a,/3, Ao, Ml, M 2 ) > 0 and r 2 = r 2 (d, ce,/3, Aq, Mi, M 2 ) G (0,ri] such that for every 

zq G dD and r G (0,r2], 

dxd9b,r{^) ^ forxeDn B{zo,rl2). 

OD\X) 

Proof. Without loss of generality we assume zq = 0. Let r 2 G (0,ri] to be specified later. For 
r G (0,r2], fix X G iA n H(0,r/2). By (15.ip . xq := p{x)/2y/l + Aq < 5d{x)/2 < r/4 < r2/4. Set 


X = X + 2p{x)ed and x = x — 2p{x)ed. 

Observe that i?(x, xq), H(x, xq) C D~^ and B{x,rQ) C D~{0,+oo,2r). By (I4.30p and (|4.3p . we have 

dxa9b,r{x) = /_ gb,r{z)dx^KLx,To)i^^^)dz 

JB{x,ro) 


> 


f B(x,ro) 


, 9b, 


dxaGBix,ro) y)J^{y, z)dy 


' B{x,ro) 


9b,ri^) 


B{x,ro) 

(/ 

\J B{x,ro)xB{x,ro) 


dz 


dxaGB{x,To){x,y)SlG'^B{x,ro)iy^w) J\w,z)dydw 1 dz. 


(5.4) 
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Let A := l/ro and Bi := Ai?(x,ro). By scaling property, Lemma 14.41 and Lemma 13.21 we have 


L 


B(x,ro) JB{x,ro) 


dxaGB{x,ro) {x, y)*S’^G'^(,r,r-o) (?/> w) j\w, z)dy dw 


< 


A'^+i [ [ 

Qx^Gbi Ax, y)Sy>-G\ {y, w) 

JBi JBi 


c^^d+i-a+gJ 

^ GbiAx,w)A^{w, \z)dw 
Bi 

^iA1-o+/3 f 

GB(x,ro){x,v)jAv,z)dv 

Jb( 

^,ro) 

2ciA^-“+^ [ 

GBix,ro)A^Aj\v,z)dv 

Je 

{x,rQ) 

2ciro^+“ ^K^B{x,ro)A^A- 


Here ci = ci(d, a, 13, Mi, M 2 ) > 0. Thus we can continue the estimate in (|5.4I) to get 
dx^9b,r{x) 


> 


/_ jhA^) 

J B{x,ro) A 


B{x,ro) 


dxaGB{x,ro){x,y)j\y,z)dy dz 


2ci 0-/3 


-n 

ro 


' B{x,ro) 


9bA^)KB{x,ro)A^Adz 


/_ ,9bA^) 

J B{x,ro) A 


2ci 0-/3 


B{x,ro) 


dxaGB{x,ro)ix,y)J {y,z)dy]dz - A 9b,r(x). (5.5) 


ro 


Note that by 

dxaGB{x,ro){x,y) 

= 2^-“7r-'^/2p p f y£jlA£r-?(r? - U, - 


\y-x 

I / 7 \yd ~ ^ / \ 

+ (a — 0 ')j^^——^GB{x,ro)\^^ y)' 


2^0 (^0 -\y- {roAo -\y- x\A + \y - x\^) 




(5.6) 


Obviously dx^GB{x,ro){x,y) is anti-symmetric in y with respect to the hyperplane H := {y : 
Vd = Xd}. For every z € B{x,ro)\ define hx{z) := fB(x,ro) 9xaGB{x,ro)ix,y)A{y,z)dy- By ([53]) we 

- C 

have for a.e. z € B{x,rQ) , 

hx{z) = / dx^GB{x,rQ){x,y) [A{y,z) - A{y,z)\ dy 

J {y^B{x,ro),yd>Xd} ^ 

= I dxaGB{x,ro){x,y) (A{y,\z-y\) - f{y,\z-y\)] dy, (5.7) 

J {y€B(x,ro),yd>xd} ^ ' 

where y = {y,2xd — yd)- We observe that the right hand side of (|5.7I) is antisymmetric in 2 : with 

respect to the hyperplane %. Recall that A{y,r) is non-increasing in r. Following from this, the 

monotonicity of 9b,r{x) in Xd and the Harnack inequality, we have 

/_ gb,r{z)hx{z)dz > / gb^r{z)hx{z)dz = / gbAz)hx{z)dz 

JB{x,rA J B{x,ro)UB{x,ro) J B(x,ro) 
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hx{z)dz 


(5.8) 


> C2gb,r{x) I 
Jb\ 


' B{x,ro) 

for some C 2 = C 2 {d,a, P, Xq, Mi, M 2 ) > 0. Note that by (15. 6 p and ()5.7I) . for a.e. z £ B{x,ro), 


hx{z) > {d-a) 

= {d — a) 


{y£B{x,ro)-.yd>Xd} 


Vd-Xd , 

\y - x\ 


GB{x,ro){x,y) [j {y, \z - y\) - j {y, \z - y\)] dy 


I B{x,ro) 


^y- J| 2 GB{x,ro)ix,y)j^iy, jz - yl)dy > 0. 


Therefore by the scaling property of the Green function GB[x,ro)^ 
hx{z)dz 


J B(x,ro) 
> {d — a) 

d — a I 


B{x,ro) JB(x,ro) 


{x, y)j'^{y, z)dy dz 


> 


M- 

d — a 


2 r-Q |u; - X 

I I I I 




[ [ — ^^4^Gb(o, i)i0, w - x/ro)J{\v - w\)dwdv 

J\v-ZL \<1 J\w-ZL\<i \w - x/ro\^ 


roM 2 7|„_J.|<i y|^_jiL|<i \w — X 

I ro I I ro I 

d — a 


/ / ■^Gb(o,i){0,w)J{\v- w + -—-\)dwdv 

J\v\<l J\w‘ “ 


roM2 J\x\<i J\w\<i kP 

d-a f f Wd 
roM2 J\v\<l J\w\<l 


ro 


w\ 


GB{o,i)i^^ru)J{\v - W + A 1 J 1 + Xled\)dwdv =: (5.9) 


with C 3 = co{d, a, M 2 ) > 0. It follows from ()5.5p . (15.6p . (j5.8p and p5.9p that 

( C 2 C 3 - 2cirn I gh.r(x) > — ^ 


dxdgb,r{x) > ^ (c2C3 - 2ciro gbAx) > (c 2C3 “ 2ci (r2/4)" gbAA- (5-10) 

The lemma now follows from (15.101) by setting r 2 so small that 2ci (r 2 / 4 )“~^ <£ 203 / 2 . □ 


Lemma 5.4. Let r 2 G (0, ri] C (0,1] he the constant in Lemma 1 5. ,11 There is a positive con¬ 
stant Cij = Cn{d,a, I3 ,Xo,Mi,M 2) such that for every r G ( 0 , 02 ], there is a constant 03 = 
£ 3 ( 4 , a,/3, Ao, Ml, M 2 , r) G (0, r/2) so that for every zq G dD and every non-negative function 
f that is regular harmonic in D Pi B{zo, 2r) with respect to and vanishes in n B{zo, 2r), 


\dx^fix)\ > Gn f^^\ for X e DriB{zo,r 3 ). 

5d{x) 

Proof. Without loss of generality we assume zo = 0. For r G ( 0 , 02 ], hx an arbitrary x G H n 
i?(0, r/(2Y^l + Aq)). Let Zx G dD be such that \x — Zx\ = p{x). Define c := YrmB^y^za, f id) / 9b,r{y) 
and u{y) := cgb^rid)- Obviously B{zx,2>rA) C B{fd,2r), and thus /, u are harmonic in D n 
B{zx,3r/2) andyanishm D^r\B{zx, 3r/2). Since Ad)/fid) = /(y)/'“(y) = 1) 

by Lemma IHTTl for any y G D P B{zx,3r/A), 


u{y) 

fiy) 


V 


fid) 

Ad) 


< Cl 


\y - Zx 

r 


71 

< C2 


(5.11) 
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for some positive constants Cj = Ci{d,a, P, Xq, Mi, M 2 ), z = 1,2. Consequently, 

(1 + C 2 )"V(y) < u{y) < (1 + C 2 )f{y) for y € Dn B{zx, 3r/4). 


(5.12) 


Note that x € D D B{zx,‘X>r/^) since p{x) < 5d{x)\/1 + Xq < r/2. By Lemma [531 and (15.121) . we 
have 


dxj{x) > dx^u{x) - \dx^{f - u)(x)| > C 3 

> - |9„^(/-m)(x)| 

dD(x) 


u{x) 

' 6 d{x) 


- \dxM-u)ix)\ 


(5.13) 


We assume p{x) < 3r/128. Set v{y) := f{y) — u{y) and ^ := 2p{x). Let p G (16p(x), 3r/8) to be 
specified later. In the rest of this proof, we set Di := D~^{zx,^,^) and D 2 := D^{zx,p,p)- Then 
B{x,5d{x)) C Hi C II 2 C D B{zx,‘2>r/ a) and 6d{x) = 5dx{x)- Define V{y) := lE^ |u|(X(!^^) . 
Clearly V is regular harmonic in Hi with respect to and |u(y)| < V(y) for all y € By 
Theorem 01 we have 


\dxAx)\ < \dxAA\ + \dx^{V - z;)(x)| < C 5 


doAx) Sd{x)' 


(5.14) 


We aim to estimate V{x). Note that 
V{x) < eJ|z;|(W(!„ ):W^„ GH 2 


■ E. 


^ € D 


^2 +E. 


^ e H 


= : I{x)+ II{x)+ III{x). 

By (15.lip , for any y ^ D 2 C D D B{zx, 2p) C H n B{zx, 3r/4), we have 


\v{y)\ = f[y) 


u{y) 


f{y) 


- 1 


< ci/(y) 


\y - Zx 


71 


< Cg 


p\'n 


fiy)- 


Thus 


Ax) < eg E,^, = cg yyy fA. 


(5.15) 


Let Ax D he such that p{Ax) = \Ax — Zx\ = p/lh. Define H 3 := B{Ax,p/l(}i/l + Xq). We 
observe that H 3 C H ft B{zx,p/2) C D 2 and Hi C H n B{zx,p/A). For any y G H 2 D supp/ and 
z G Hi, we have \y — Ax\ > 7pjlG, \Ax — z\ < 5p/16, and 


\y - z\>\y- Ax\ - \Ax - z\> -\y - Ax\. 


(5.16) 


If we let A := l/diam(Hi), then 
(j2.5l) and Lemma ESI we have 


= diam(Hi)" ^||5||oo < {^Pix))^^ ^Mi < Mi. Thus by 


G\)Ax,z) = X‘^ “G^^(Ax, Az) < cylx — 2 :|" z G Hi (5-17) 

for some constant C 7 = C 7 (d, a, j3. Mi) > 0. So by (j5.16p and (I5.17|) . for any y G Hf fl supp/, 

KAA^,y) = [ G],Ax,z)j\z,y)dz<cjM 2 I \x - zA-A^A\y - z\)dz 

Jdi Jd^ 
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< 

r\-i 


C7M2J^°i\y-A,\) 


\x — z 


la—d 


dz<C 7 M 2 rJ"°(ly-A,l). (5.18) 


On the other hand for any y € ^ supp/ and z G II 3 , we have \y — z\ < \y — Ax \ + \Ax — z\ < 

12\y — Ax\ll. Thus by Lemma [32] and (II. 6 p 




ID 3 


> M, 


GD3{^x,z)J^{z,y)dz 


X A'q\^r°{\y - A^\). 
Combining (I5.18P and (I5.19p . we have 

Consequently, by (15.201) , Lemma 15.11 and (15.ip , we have 

II{x) = f f{y)K’[)^{x,y)dy f 

JDinSUBDf V Jl 


~ ^ ^D3 i^x, y) for y G ZI 2 n supp/. 


f{y)KDi{x,y)dy < — / f{y)Kj^^{A^,y)dy 

1 JD^nsuppf 

< ^ I f{y)K^D3{A.,y)dy = ^f{A^) 


'D^nsupp/ 

r 


r Jd- 

~ 7,” p(l)“-T5 ■'' ' 


7]^ 


p{x, 


72 


t 72 


/7) 


Similarly we can prove that 




■ 1 I 2 


Combining p5.15p . (I5.2ip and (I5.22p . we have 

V{x) < (^(3 (Jj-Y + C8 
Thus by (I5.13p . (I5.14p and (|5.23l) we have 


p{x] 


72 


fix)- 


dxjix)> [c4-C5(ce— + cs 


yT'l p(x)T'2\\ fix) 


6d{x) ' 


(5.19) 


(5.20) 


(5.21) 


(5.22) 


(5.23) 


(5.24) 


Let rj = lQpix)'^'^/Gi+'y 2 ) ^ lemma now follows from (15.241) and (15.ip provided we choose r-^ small 
enough such that C 5 (cglO'^^r”'^! + cglO”'*'^) (r 3 (l + < C 4 / 2 . □ 

Proof of Theorem II.5t The upper bound in (|1.8p was established more generally in Theorem 
11.41 The lower bound follows from Lemma 15.41 and the inequality |V/| > \dx^f\- □ 


Remark 5.5. Taking b{x, z) = e € (0, M 2 ] in Assumption 1, we get from Theorem ll.4l and Theorem 
01 the uniform gradient estimate for mixed stable processes, which in particular recovers a main 
result of [22] on gradient estimates. 
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6 Examples 

In this section, we give some concrete examples where Assumptions 1 and 3 hold. 


Example 6.1. If b{x, z) = l{| 2 |<ci} for some ci > 0, the jumping kernel of the corresponding Feller 
process is 


J^{x,y) = 


A{d, —a) ^ A{d, —f5) 


\x — y\d-+^ ' |j; _ y|d+/9 

In this case is the independent sum of a symmetric a-stable process and a truncated symmetric 
/3-stable process, and Assumptions 1 and 3 hold with Sq = 0 and ^'(r) = I{^<c^}, respectively. 

More generally, suppose b{x,z) = hi{x, z)l{\z\<ci} for some ci > 0 and a bounded function 
6 i(x,z) on X that is symmetric in z and is bounded between two positive constants. Then 
Assumption 1 holds with eg = 0- D 


Example 6.2. If b{x,z) = 1 + \z\^ fo'^ some C 2 > 0 and 0 < 7 < /3, the jumping 

kernel of the corresponding Feller process X^ is 


J^{x,y) 


■^{d, -a) ^ A{d,-/3) 

\x — y\d+a _ y^d+y 


M((i,- 7 ) ^ 
\x - y\‘^+'r 


In this case is the independent sum of a mixed-stable process and a truncated symmetric 7 - 
stable process, and Assumptions 1 and 3 hold with eq = 1 and ip(r) = 1 -|- r^~'>'I|y.<e^j., 

respectively. 

More generally, suppose b{x, z) is a bounded function on x that is symmetric in z and is 
bounded between two positive constants. Then Assumption 1 holds with eg = 1. □ 

Example 6.3. We consider the following stochastic differential equation on 


dXt = dYt + C{Xt-)dZt, 


( 6 . 1 ) 


where T is a symmetric a-stable process, Z is an independent ;0-stable process with 0 < /3 < a, and 
C is a bounded Lipschitz function on M'^. Using Picard’s iteration method, one can show that for 
every x € M'’*, SDE (16.11) has a unique strong solution with Ag = x. The collection of the solutions 
{Xt,¥x,x G forms a strong Markov process A on Using Ito’s formula, one concludes that 
the infinitesimal generator of A is with b{x,z) = \C{x)\^. If there exists cg > 0 such that 
|C(x)| > C3 for X G M'^, then our Assumption 1 holds with eg = 1. □ 
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